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Kination dominated quintessence models of dark energy have the intriguing feature that the relic abundance 
of thermal cold dark matter can be significantly enhanced compared to the predictions from standard cosmology. 
Previous treatments of such models do not include a realistic embedding of inflationary initial conditions. We 
remedy this situation by constructing a viable inflationary model in which the inflaton and quintessence field 
are the same scalar degree of freedom. Kination domination is achieved after inflation through a strong push or 
"kick" of the inflaton, and sufficient reheating can be achieved depending on model parameters. This allows us 
to explore both model-dependent and model-independent cosmological predictions of this scenario. We find that 
measurements of the B-mode CMB polarization can rule out this class of scenarios almost model independently. 
We also discuss other experimentally accessible signatures for this class of models. 



I. INTRODUCTION 

The discovery that the universe is dominated by dark energy strongly suggests that the predictions of standard cosmology 
should be reevaluated. An intriguing possible explanation of the nature of dark energy arises within the quintessence paradigm, 
in which the dark energy takes the form of a slowly evolving scalar field (see e.g. d 0. [H 0, H[|). This scenario is extremely 
difficult (if not impossible) to test directly in collider experiments, since quintessence models generically require the quintessence 
field to have gravitationally suppressed interactions with the fields of the Standard Model (SM). 

However, if the dark energy is in the form of quintessence, the presence of the quintessence field can modify the cosmological 
evolution and lead to significant departures from standard cosmology. A striking example is the possible interconnection of dark 
matter and dark energy within quintessence scenarios la 13 l2l llOIl • As first pointed out by Salati [6], the freeze-out of thermal 
relics can be strongly enhanced in scenarios in which the energy density is dominated by the kinetic energy of the quintessence 
field (kination domination) during the time of freeze-out, but dilutes away by the time of big bang nucleosynthesis (BBN). 
(Related scenarios were also suggested before by UHO-) Such kination dominated freeze-out scenarios are then consistent 
with standard cosmology and predict that the standard relic abundance computed from the parameters extracted from collider 
measurements will be mismatched from the relic abundance deduced by observational cosmology. This has implications for 
TeV physics models with thermal dark matter candidates (e.g. models with low energy supersymmetry such as the minimal 
supersymmetric extension of the SM (MSSM), technicolor models, models with large/warped extra dimensions, or certain 
classes of little Higgs models), which will be probed at the LHC and other experiments in the foreseeable future. 

In most of the previous discussions of this class of scenarios H 0,111 13], the initial condition that the quintessence field kinetic 
energy density is the dominant component was put in as an ansatz, without a complete picture of the inflationary dynamics. In this 
paper, we address this issue by constructing an inflationary scenario that dynamically leads to a kination dominated quintessence 
period. This yields robust predictions which can be used to experimentally support or rule out this class of scenarios. 

To construct viable inflationary models which lead to kination domination, the following constraints must be satisfied: 

1 . The energy density in a coherent quintessence field must dominate over radiation after the end of inflation. 

2. The quintessence field must be kinetic energy dominated. 

3. The inflaton potential must satisfy the usual requirements of a sufficient number of efoldings, the right amplitude of density 
perturbations, and a nearly scale invariant spectral index (with a slight preference for a red spectrum [ 13J). 
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4. There must not be too much reheating in the phase transition at the end of inflation (when the quasi-de Sitter phase ends), 
such that the ratio of the kination energy density to the radiation energy density can be large. 

To satisfy these constraints, we build a model in which the quintessence field is the infiaton field $, which has sufficiently weak 
couplings to the SM fields at the end of inflation such that most of the energy density responsible for inflation gets converted to 
$ coherent kinetic energy in a runaway potential. Radiation domination is achieved because the coherent kinetic energy dilutes 
as 1/a 6 , while the suppressed radiation produced at the end of inflation gets relatively amplified since it dilutes as 1/a 4 . The 
choice of the infiaton as the quintessence field in kination domination is natural in this framework, since the infiaton possesses 
the required qualities for the kination domination construction: energy dominance and coherence. Other quintessence models 
have been constructed in which the infiaton is the same field as the quintessence field (see e.g. 11141 Il5l llq, Il7l ll8IO . but we 
treat the reheating more precisely and derive new predictions in the context of kination-dominated quintessence models relevant 
for dark matter. The work tlOfl considers the connection between dark energy and dark matter, but it differs from the present 
treatment in that their scenario is model specific and has multiple inflationary phases. 

In addition to the prediction that the relic abundance inferred from collider measurements can measurably disagree with the 
cosmologically inferred dark matter relic abundance (in the context of a thermal freeze-out scenario), many experimentally ac- 
cessible cosmological predictions can be made by embedding the kination domination scenario in inflationary models. The most 
model independent nontrivial prediction is the absence of a measurable tensor perturbation induced B-mode CMB polarization. 
Hence, the next generation of CMB experiments can falsify this class of models. 

Other predictions include a shift in the peak of the gravity wave signal from the electroweak phase transition, shifts in the 
verifiable leptogenesis/baryogenesis scenarios whose out of equilibrium ingredient is furnished by H, and shifts in any residual 
annihilation effects. In general, the indirect dark matter detection signals are enhanced in this scenario 11911 . In particular, it 
would be interesting to see whether one can explain the positron excess as observed by HEAT and other experiments near 7 GeV 
within this scenario 1 20], since most attempts to explain this excess in terms of dark matter annihilations require sufficiently 
large cross sections that the relic abundance would be negligible with standard cosmological assumptions. 

The order of presentation will be as follows. We begin by analyzing both analytically and numerically the constraints on 
quintessence potentials to achieve a period of kination domination, and look for equation of state signatures. In SectionUm we 
construct a class of inflationary models embedding the kination domination scenario and discuss a robust prediction that can 
observationally falsify this class of models. In SectionHVl we discuss other observationally accessible cosmological predictions 
which may corroborate this class of scenarios. We then summarize and conclude. In the Appendix, we provide the details of the 
particle production computation for the unusual reheating scenario associated with this class of inflationary scenarios. 

Throughout this paper, we use the convention of M p to denote the reduced Planck mass of approximately 2.4 x 10 18 GeV. 



The standard procedure in constructing scenarios of quintessence dynamics is to focus on potentials, which is a sensible 
approach since a negative equation of state requires potential energy domination. However, when considering quintessence 
dynamics with a period of kination domination, it is more natural to focus on the behavior of the field velocity, as this is the 
quantity which characterizes the energy density of the quintessence field. Here we develop a formalism to map the desired 
behavior of the kinetic energy to a class of scalar potentials. Our results show that the allowed quintessence potentials are 
not severely restricted by the requirement of a period of kination domination, since the time at which the equation of state is 
close to —1 is necessarily much later than the time of dark matter freeze-out due to the strong constraints from BBN. In what 
follows, we will not restrict ourselves to "tracker" models, in order to separate the difficulties of constructing good trackers from 
the constraints imposed by kination domination. It is worth noting that kination dominated initial conditions can be obtained in 
tracking potentials, as discussed in [ 8 , 2 1 ] , although clearly there will be constraints depending on the exact form of the potential. 

To carry out the engineering of quintessence potentials which match the desired kinetic term histories, we begin with the 
familiar equations of motion for the quintessence field $ and the Friedmann-Robertson- Walker (FRW) scale factor H = a/a: 



in which prm = Pr + Pm corresponds to the energy densities of radiation and matter. Using the following definitions, 



II. MAPPING KINETIC BEHAVIOR TO QUINTESSENCE POTENTIALS 



$ + 3i?$ + V r/ ($) = 



(1) 




(2) 



= ma 3 ($) 



(3) 



/(*) = Ht) 



(4) 
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(5) 



we rewrite Eqs. (03 and (|2|i as follows: 



/ 2 ($) /(*). 



(6) 



= h?l^p(<S>) + V(<P)+p RM ( q (<S>)). (7) 

In the above, we have assumed that $ is a single valued function of $ (which would exclude e.g. the situation of oscillations). 
For the scenarios considered here this is typically not a good assumption throughout the entire evolution, since the kinetic 
energy is usually large enough that $ can overshoot the minimum and eventually hit the "other side" of the potential before the 
quintessence equation of state reaches — 1. However, the condition on the quintessence potential during kination domination 
does not change even if this behavior is properly accounted for since the "bounce" occurs long after the kination period is over. 

We can use Eqs. (|6]l-(|7]i to solve for q(Q) and /(<£>) if V (<&) is given; alternatively, we can solve for V($) and q(&) if the field 
velocity function /(<£>) (i.e. the kinetic energy) is specified. For generic potentials, there is no obvious obstruction to achieving 
a period of kination domination which leads to a period of potential energy domination, as /($) can be chosen to vanish as $ 
approaches a particular asymptotic value, and the Hubble friction naturally allows for /($) to vanish without V also vanishing. 

In this paper, we restrict our attention to the specific case of the quintessence field velocity function with the initial behavior 
(during the time relevant for dark matter freeze-out) of 

/«/ie«-*oci, (8) 

in which fi and qi denote initial values of the functions /(<£>) and q(Q). This corresponds to the kination regime in which the 
universe is driven by quintessence kinetic energy, with V' / f 2 playing a subdominant role compared to /'// in the equation of 
motion for $ (i.e., the force term is subdominant to the Hubble friction term). The reason of course for this restriction is that 
kination domination gives rise to an energy density which dilutes as 1/ a 6 , which allows the quintessence energy to be important 
during DM freeze-out but disappear by the time of BBN as required by phenomenology. 

Let us now study what a potential that leads to the kination behavior of Eq. (0 looks like by constructing it using Eqs. © and 
(IT). To obtain a closed form solution, we take the ansatz of neglecting V'// 2 , in which case we obtain 

<Z = % -lnffl (9) 



^/ 2 ($)-p flM (9($)) + ^[/'($)] 2 . 



(10) 



One possible procedure is then as follows: choose any desired behavior of / and obtain using Eqs. (0 and ( fTOb a corresponding 
expansion history (specified by q) and the potential V, then check that V / f 2 can be neglected in the equation of motion. The 
validity of neglecting V / f 2 can be rewritten as 



-1 



1 dpRM 
P dq 



97 / 



< 1. 



(11) 



One trivial way to make V / f 2 negligible is to have V much smaller than any other energy component and to have V be a 
smooth function. Another way to satisfy Eq. ( fTTb is to set its left hand side equal to an arbitrary small function h(Q) <C 1 and 
solve for /(<&). For example, if the universe is radiation dominated, we can write 



Prm (q) = Pi 



a,i \ 4 
a 



Pie 3 



to obtain the following equation: 



/"($) - 67 



n ) 



1/3 



Kfi 



9-) 



(1 + &(*))/($) = a, 



(12) 



(13) 
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in which 

K = Pi /f?K^^, (14) 

with 77$ = p$//0 7 |t=imcV^ an d fi = /(^i) ~ 6 x 10 3 y / 77iTGeV 2 (the initial field velocity) as defined at a temperature of 
approximately 1 GeV. Note that since h(Q) <C 1 is arbitrary, nearly any smoothly varying f(Q) can be obtained, which in 
turn means any smoothly varying shape for the potential can be obtained even with the constraint of Eq. (0. (This corresponds 
of course to the case in which the V energy density is negligible during kination domination.) To obtain intuition, if we set 
h($>) = h = constant <C 1 and define y = /($) / fi and x = we can rewrite Eq. ( fT3l l in terms of a first order equation: 



^ = -yW/3 + |(l + ^ ) tf 2 + C J (15) 

in which C is an integration constant. Note that we have chosen the negative squareroot to coincide with the convention in this 
paper that $ is moving in the positive direction during kination domination (recall y is decreasing by construction). 

Although Eq. (fTBI l can be used to write the solution as a single integral, the phase space of solutions can also be visualized by 
examining this equation. As long as C > 0, y moves towards y = (y ~ 2 x 10 -10 at the time of BBN). The existence of a 
large class of solutions given by C > is not surprising given that the Hubble expansion generically provides friction for the 
field velocity /. The time it takes for the field to achieve a final velocity // can be written as an integral: 



dy 
y dx 



As a consistency check, it can be seen from Eq. (fTBI l that generically it takes an infinitely long time to achieve // = (an infinite 
time is necessary since we are engineering / oc a~ 3 ). Although a larger C seems to correspond to a shorter time to achieve the 
desired //, it also implies that the potential is larger, which places severe constraints on C. More explicitly, it can be shown that 
the potential of Eq. ( fTOb implied by Eq. ( fTBI l is of the form 

where C is then explicitly seen to control the cosmological constant. (If C < 0, the cosmological constant will be negative, 
which will cause the universe to eventually contract. This is consistent with the C > condition described above.) Note that 
the potential is independent of k except through y(x). It can be seen from this expression that if h($>) were not a constant ho, a 
richer shape of potentials can easily be attained (i.e., it is easy to show that the y dependent term generalizes to ff J dyyh when 
h is not a constant). Furthermore, we see the condition that V($>) be a subdominant energy component to / 2 /2 implies that 
C <C 1. Specifically, if Cff/9 is identified as the cosmological constant that persists to today, one would impose the constraint 
C < 9flAPc/ ff ~ 6 x 10 -54 /?7$. If we insist on this form of the potential until BBN, we have a bound of 

10 -i8 

C< (18) 

even if this piece of the potential were piecewise glued to other functional forms for the potential once the evolution has persisted 
past the kination dominated period. Since C is generally required to be small during kination domination, we can now solve 
Eq. (fTBT ) exactly neglecting C, which leads to the expression 

y(x) - + 1 fe )3/2 exp(-3v/l + h Q (x - C 2 )/V2) (l - 8k cxp[V2(l + h )(x - C 2 )]) 3 , (19) 

in which C2 is another integration constant that is specified by y(xi) = 1: 



n V2 1 

C 2 =Xi+ In 

VT+ho 



2(Vl + ho + \/l + h + 2k) 



(20) 



As expected, initially the radiation energy density encoded in n is unimportant, but once the radiation catches up with the 
quintessence kinetic energy, the potential must compensate to (artificially) maintain the 1/a 6 behavior for the quintessence 
energy density while keeping the left hand side of Eq. ( fTTT i. h(Q), a constant. Once again, it is important to note that the shape 
of the potential dictated by Eq. ( fT9l is not fundamental because this solution is only valid for a constant h(&). 
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The main lesson from the discussion thus far is that potentials can be chosen to not interfere with 1/a 6 behavior of the 
quintessence energy density even past the point in which the radiation starts to dominate. (As discussed earlier, this follows 
from the fact that the solution to the $ equation of motion with V'(Q>) = is $ oc a~ 3 .) To obtain a potential that is manifestly 
independent of k, we can choose implicitly a non-constant h($>). For example, inspired by Eq. ( fl"9l ), we can choose the potential 
to be 

V = fi A Pc[l + 6cosh 2 (A$)], (21) 

in which 



A ^VT (1 + 6) ' (22) 

with b <C 1. This form of the potential results in V'($>) being unimportant during the $ evolution until the temperature reaches 

T ~ 10- 10 6 2 / 9 (1 + b) 1 ' 9 exp (-10VIT& + t,- GeV, (23) 

where the initial value of $ was parametrized as — n/A. The values n = 30, b = 10~ 6 , and 77$ = 1 result in a temperature 
during matter domination but close to the matter-radiation equality (i.e., long after the end of BBN). Therefore, this potential 
provides the desired kination behavior. We will be using this form of the potential for the rest of the paper. 

To check the stability of the background scalar field solutions, we set $ = s (t) + a~~ 3 / 2 5(j)(t,x), in which <^s(t) is the 
background solution, and note that upon neglecting the metric fluctuations to a leading order approximation, the equations of 
motion for the field fluctuations take the form: 

64> - ^(d 1 64>) 2 + [V"(ct>B) + ^P]S4> = 0, (24) 

a z 4 

in which P is the pressure of an ideal fluid (our stress tensor approximation). The condition for stability on all length scales is 

V"{4> B ) + \lP > 0. (25) 

Therefore, if V"(<j>B) > (which is true for example in the case of constant /i(3>) in Eq. ( fTTI i if ho > or in the case of Eq. (12 li t 
if b > 0), the solution will always be stable for the case of P > relevant for kination domination scenarios. 

The viability of the quintessence picture also requires that the background solution tends toward a potential energy dominated 
regime with V($) ~ QaPo m which p c is the critical energy density today. Since the equation of state can be written as 

l-2V/p + l PR /p 
W l + 2V/p + 2(p M + PR )/p> {0> 

it is straightforward to see that as long as V asymptotically dominates the energy density as the universe expands, the kination 
scenario naturally leads to the desired late time quintessence behavior of w — > —1. The phenomenological requirement is that 
w w —1 by redshift of about z ~ 1. 

A final issue to address is the underlying dynamics which leads to the required initial conditions to achieve the desired /j. We 
will return to this question later and will discuss interesting observational consequences, but first analyze an explicit example 
numerically to illustrate the possible resulting equation of state. 

In solving the quintessential cosmology specified by Eq. (BIT ) numerically, we choose the following representative parameters: 
n A = 0.72, fl M = 0.28 (including baryons), fl R = 4.6 x 10~ 5 , and H = 73 km/s/Mpc. We further set b = 10~ 6 , 
A<i>(ttinitiai) = {—30, —20, — 10}, r7$ = p$/p 7 \t=imcV = 0.5, and A<j>(ii n i t i a i) > (the magnitude is fixed by 77$). We also take 
the effective degrees of freedom to evolve approximately as 

{90 T > 1 GeV 

60 1 GeV > T > 0.1 GeV 
10.75 0.1 GeV > T > 10~ 4 GeV (27) 
3.36 10~ 4 GeV > T. 

Although the relation g*(T) w g*s(T) breaks down at late times, these subtleties do not affect the main results for kination 
domination. In evolving the coupled differential equations for {a(t), $(t)} as given in Eq. (Q]) and Eq. ©, we take 

f*\ i, a ^ g*(W V /3 ^(today) V 
PR (t) = Pr (today) — — y — -— (28) 

VS*(T'approx(t))y V a ( f ) / 
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log(T/GeV) 



Figure 1: <J> as a function of T. The evolution for A$(i t initiai) = { — 10, —20, —30} is depicted by the curves denoted by Models A, B, and C, 
respectively. The little jump in temperature at T w 10~ GeV corresponds to the change in g, at that temperature; a much smaller jump can 
be seen at T « 1 GeV. 




Figure 2: The dark energy equation of state w as a function of redshift z for the same cases shown in Fig. [TJ 



and 



PM{t) 



,'a(today)Y 
PM (today)(A_ij 



(29) 



This neglects the fact that pm has a nontrivial time behavior due to annihilations; however, the annihilation corrections are 
negligible by the time pu becomes a significant component of the energy density. We are also using an approximate temperature, 



^approx {tj — ^today 



^today 



(30) 



in 3* in Eq. ( |28l l. This should be sufficiently accurate at the current level of approximation, since is a very flat function except 
in critical regions of Eq. d27l i where the transition temperatures are only accurate to an order of magnitude. 

The resulting evolution of $ is shown in Fig.Q]as a function of temperature. For A$(ttmitiai) = —30, which we denote as Model 
C, the field does not reach the origin (the minimum of the potential) by today. However, when the initial value of A$ is less 
negative, the field overshoots the minimum and climbs up the potential for $ > 0, and eventually falls back towards the origin. 
An inspection of Fig. [2] which shows the dark energy equation of state as a function of redshift, demonstrates that acceptable 
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Figure 3: The equation of state w of the $ component as a function of redshift z for large redshifts for Model A (A<&(ttniitiai) — —10). 

phenomenology can be obtained for the "overshooting" case of A<f>(i t i n j t i a i) = —20 (Model B), but not for A<f> (initial) = —10 
(Model A), where there is still too much kinetic energy by redshifts of less than one. In fact, the equation of state w for Model 
A actually reaches —1 during the "bounce" and then a time period in which w = is sustained until the kinetic energy is 
finally dissipated. This nontrivial evolution of the equation of state for larger redshifts than those shown in Fig. [2] can be seen 
in Fig. [3] Furthermore, it can also be shown that the scalar field energy density always dominates in Model A, which makes 
such scenarios incompatible with the successes of CDM large scale structure phenomenology. Clearly, even with exactly same 
potential function V(Q) and initial value of quintessence kinetic energy, the dark energy phenomenology crucially depends on 
the initial value of $. 

In all of these cases, the initial condition of | A<I> (initial )| > 10 corresponds to situations in which the initial displacement is 
trans-Planckian. This is not attractive from an effective field theory point of view since the theory then is sensitive to all powers 
of the field operator; however, such initial field values are typical in quintessence scenarios B22I1 . Indeed, although quintessence 
is genetically not attractive as an effective field theory candidate, its classical dynamics may still be useful to parameterize 
vacuum dynamics, since effective field theory has grossly failed as far as understanding the vacuum structure is concerned, i.e., 
in giving a plausible explanation for the cosmological constant problem. In fact, we have implicitly been assuming throughout 
this paper that the cosmological constant problem has been solved by some unknown mechanism which leaves the quintessence 
dynamics responsible only for the nontrivial vacuum energy dynamics. 

The main lesson to be learned from this exercise is that in viable kination scenarios, the dark energy equation of state can 
exhibit a wide range of behavior both because of the shape of the potential and because of the unknown initial conditions. How- 
ever, a typical behavior more likely in quintessence scenarios with a period of kination domination than in generic quintessence 
models is the "bounce" or turnaround behavior, in which $ bounces off the "other" end of the potential barrier and slowly rolls 
toward the minimum. With one bounce, the behavior of the Model A case in Fig. [2] in which dw/dz < today is typical in 
viable models. The complement, dw/dz > 0, is less generic for large initial kinetic energies, but may be possible in models 
in which the kinetic energy during the kination period is sufficiently small and tuned appropriately. Multiple bounce scenarios 
may also be possible in kination scenarios, but they require a special conspiracy between the shape of the potential and initial 
conditions, as exemplified by the unsuccessful case shown in Fig. [2] 

In the present model, we can compute the minimum initial displacement | A $ (^initial )| to avoid the turnaround behavior. The 
total displacement of $ for the case in which V'($) plays a negligible role in the equation of motion is 

A$ w V6Mp In 

in which a(ij) is the scale factor at the time of the initial position of $, and the second term on the right hand side of Eq. (|3TT i 
arises from the movement of $ after BBN. For typical situations in which a(i BBN )/a(ij) = a(t BBN )/a(tp) > 10 3 , where a(tp) 
is the scale factor at the time of freeze-out, we find 

A$ > [17 + 3v^i]A/ p . (32) 

This result implies that most of the displacement of $ occurs during the kination period and not after BBN, and that the numerical 
result of Fig.[T]is reasonable. In addition, the field displacement is generically transPlanckian (as discussed previously). 




o(£bbn) 
a(ti) 



(31) 
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Having discussed the consequences of having large quintessence kinetic energy for late time cosmology, we now consider 
how such initial conditions can be established within the context of inflationary cosmology. 



III. INFLATIONARY COSMOLOGY 



In this section, we address the issue of how kination domination might be achieved within theories of cosmological initial 
conditions such as inflation. To demonstrate an existence proof of a viable scenario, we construct an inflationary model in which 
a single scalar field plays the role of both the inflaton and the quintessence field. We will not only show that kination domination 
can be achieved after inflation, but that sufficient reheating can be achieved as well. The main benefits of embedding quintessence 
within an inflationary scenario are the resulting cosmological predictions which can be experimentally verified or falsified. We 
will see that a robust prediction of the kination scenario relevant for dark matter abundance is an absence of a measurable B- 
mode CMB polarization signal. Other predictions will include those connected to the fact that the Hubble expansion rate during 
kination domination is different from that within standard cosmological scenarios even when the radiation is in equilibrium. 

The basic idea we implement to construct this inflationary model is the assumption that the inflaton receives a kick (or a strong 
push) at the end of inflation to achieve kination domination. Indeed, since slow roll inflationary models already require a coherent 
homogeneous scalar field to dominate the energy density in the universe, the inflaton is an ideal candidate field to be converted 
into a quintessence driving kination dominated universe. By comparing a multi-scalar field system in which the inflaton energy 
density converts efficiently into a coherent scalar field kinetic energy of another field direction, it is straightforward to see that 
the only significant simplification being made by considering one scalar field is in the neglect of the acceleration of the field 
velocity vector direction in field space. 

To construct this model for kination cosmology in the context of inflation, we will assume the following degrees of freedom: 

1. Inflaton and quintessence. A real scalar field degree of freedom <i> plays the role of both the inflaton and quintessence. 
Of course, a complete supersymmetric embedding of this scenario would require a complex scalar degree of freedom. We 
neglect this detail here for simplicity, as we are first concerned with generic dynamical settings, and note that there is no 
straightforward insurmountable obstacle to extending this scenario into a fully supersymmetric framework. 

2. MSSM fields. We assume the presence of electroweak to TeV scale MSSM field degrees of freedom; the lightest neutralinos 
are stable LSP's. We will denote the neutralino LSP as \ an d other generic fields as ip. 

3. Couplings. We assume for simplicity that $ is coupled to MSSM fields only through the minimal gravitational coupling. 
Scenarios with additional couplings are of course plausible, but we do not consider them in detail here because the 
resulting constraints are highly model dependent. Furthermore, since viable phenomenology requires the quintessence 
field to be very weakly coupled to the (MS)SM fields, the minimal gravitational coupling scenario is "natural" within the 
quintessence paradigm (of course, quintessence itself has a doubtful status from the point of view of effective field theory). 

For simplicity, we will also functionally tune the potential (i.e., assume a specific form of the nonrenormalizable operators). Such 
simplifications are reasonable for this first attempt at model building of this kind, given the current incomplete understanding of 
possible UV completions of the MSSM as well as the typical difficulties of embedding inflation and quintessence in the context 
of effective field theories. Future model-building attempts will need to address this issue. 

Our scenario is predicated on the physical picture that $ receives a kick at the end of the inflationary period. To this end, we 
consider a scalar field potential with a stepfunction-like behavior, such that $ is potential energy dominated at the top of the step 
and then becomes kinetic energy dominated when $ drops off the cliff of the step. As an example, consider the ansatz 



w il APc [l + 6cosh 2 (A$)] 



Vn /31n( V ^ 



Sm, (33) 



in which (3 is a constant energy density scale controlling the slow roll properties of the period of inflation, 8 is a small constant 
inserted purely to regularize the logarithm when $ — <1> C vanishes, and is a steplike function. An example would be 

5($) = (1 — tanh[a(<£> — <£> c )])/2, which is unity for a($ — $ c ) < 1 and smoothly goes to zero after a(<& — <£> c ) > 1 for 

sufficiently large values of a. This potential is identical to Eq. (ETT) except with the addition of the Vo + f3 In ^ S(3>) 

term responsible for the inflationary period, which in turn shuts off at $ = $ c . Note that the constant S plays no important role 
but to make the logarithmic function regular, and as we will comment more explicitly later, we can choose a sufficiently small 5 
as not to change any of the inflationary analysis. Hence, it can be dropped from the anlysis of the inflationary period. 

Let us consider the constraints on the scales Vo, a, and (3 in this scenario (we will set [i = $ c without any loss of generality). 
Inflation occurs for (<E> — <£> c ) < 0, and ends when the inflaton receives a hard kick at $ sa $ c , leading to kination domination. 
During the kick at the time t e of the end of inflation, the non-adiabatic time variation of the background gravitational field 
generates particles with energy density p^(t e ) (a similar reheating scenario was considered in 112311 ). The mechanics of this is 
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explained in the appendix. If the particles that are produced have masses much smaller than the expansion rate H e ~ y/Vo/3/M p 
at the end of inflation, then 

P+ (*e) ~ %9*{T e « H e /2TT) (jA , (34) 

in which g* counts the number of light degrees of freedom. This can be seen as the situation where the all of the species which 
couple to the large $ vacuum expectation value (VEV) have decoupled and the others are lighter than H e . We have made a 
special simplifying assumption that there is a large number of light species despite the large $ VEV. (Indeed, the lack of such 
light particles can lead to a moduli problem as the VEV can induce large masses to the fields to which it couples. Note that even 
in regions where the $ VEV is zero, large finite density masses can be induced for what would otherwise be light fields ll24ll25ll .) 
Due to the kick and the sudden drop of the potential, the kinetic energy of $ just after inflation ends will be of order Vq and 
will dilute as 1/a 6 . Athough the relativistic species ip initially is out of equilibrium, the energy density can be characterized by 
an approximate temperature T e ~ (py / g*(T e )) 1 ' 4 '. As we will discuss below, tp will eventually equilibrate and the relativistic 
energy density will dilute as 1/a 4 . 

In this class of scenarios, the initial values of p^, and p$ are interconnected, such that for any given reheating energy density 
p^, there is a predicted value of 77$ = p<t>/p 7 \t=i MeV> which is phenomenologically required to be less than about 1 (at 2a) by 
the time the photon temperature is of order 1 MeV Since p$/ ' p 1 w a~ 2 oc T 2 during the kination period, 77$ is given by 

oo (1 MeV) 2 K) 

m ~ 33 ^U5 ; (35) 

P f{T e )^r7) 

Combining Eqs. (l34l and (l35l l. we find that in order to obtain a desired 77$, the inflationary energy density must be of the form 

V ~ (3.9 x 10 13 GeV) 4 % 1/2 . (36) 

A larger 77$ requires a smaller Vq, because the radiation energy density at the end of inflation is proportional to V Q 2 , and an 
increased radiation energy density corresponds to an increased scale factor growth before the temperature reaches 1 MeV. 

Eq. (l36l l is a remarkable result as it sets an approximate upper bound on Vq if a non-negligible 77$ is to be achieved. Note that 
this did not depend on the details of the inflationary model, but only on the fact that a period of kination domination occurs just 
following inflation together with reheating. Although such upper bounds have not been imposed in previous studies |6j,l8[], this 
result appears to be quite generic within a large class of inflationary models. 

An interesting ramification of this bound on Vq results from the fact that the detection of inflationary tensor perturbations in 
the foreseeable future requires Vq > (3 x 10 15 GeV) 4 (corresponding to a tensor to scalar ratio of about 10~ 4 ) [27, 28, 29, 30]. 
Hence, if tensor perturbations are detected, 77$ must satisfy the following approximate bound: 

m < 10 ^(°Y(^y>(t£l)-\ (37) 



10/ vio-v V 100 



In the above, 9 = V($jv)/V($ e ) is the ratio of the potential between the time when the largest observable scales left the 
horizon and the time when inflation ends, and r m i n is the minimum detectable tensor to scalar ratio defined as 16e evaluated 

at the 0.002MPC" 1 Hubble crossing scale, where e = -if-(V (4>)/V(4>)) 2 is the usual potential expansion slow roll parameter. 
Therefore, larger values of 77$ would be ruled out in this scenario if tensor perturbations are measured. Note that a nonzero 
77$ <C 10~ 6 causes an 0(1O 6 ?7$) change in the prediction of the relic abundance compared to the standard scenario [6, 8, 
12611 . In the event of a positive tensor perturbation measurement, 77$ would be bounded to be so small that kination dominated 
scenarios cannot be effective in changing the thermal relic abundance from the values obtained using standard cosmology (i.e., 
the connection between collider measurements and dark energy would be lost). 
The constraint for obtaining the required number of efolds is given by: 

N > In (Ve^f) + \ In ( MP**)) + m JL, (38) 
V Ho J 3 V g*s{To) J T RH 

in which Hq is the expansion rate today, and in our scenario the reheating temperature Trh ~ H e /2ir, indicating that the explicit 
H e dependence drops out of Eq. d38l ). The requisite number of efoldings is slightly larger than the usual inflationary scenario 
with reheating during coherent oscillations since there is no corresponding scale stretching that would have occurred during the 
coherent oscillations phase. As there is no explicit H e dependence, for 9 ~ O(10) the required number of efolds of inflation is 

N > 71, (39) 
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nearly independently of the inflationary potential. Note that in the absence of an unusual initial state of the vacuum, the number 
of efolds required after the onset of inflation to achieve the usual inflationary predictions is very small 13 ill . 
Since Vq is fixed by Eq. (l36l l and e is fixed by the density perturbation constraint 

F TC (0.002 Mpc- 1 ) = « 2 x 10~ 9 , (40) 



247r 2 eM4 



f3 is approximately given by 



10"'GeVi'('^)) _ (-L). (41) 



For this class of models Vb S> /3, which implies that 9 « 1 during the approximate 70 efoldings of interest since the potential 
depends only logarithmically on <f> during inflation. Nonetheless, it is still useful to keep the 9 dependence explicitly to keep track 
of the model dependence, since expressions like Eq. d37l > are not strongly dependent on the shape of the inflationary potential. 

To determine the very model dependent field value during which k « 0.002Mpc _1 leaves the horizon (corresponding to about 
68 efoldings before the end of inflation), we use the usual field integral over 1 / ^/2e($), which yields 

* - $ c « -8 x 10- 5 ( m y 1/4 (^f^) M p « $ c . (42) 

Eq. d42i > implies that unless there is fine tuning of the initial conditions, the number of efoldings will be much larger than the 
requisite number of Eq. ( |39l l. Finally, the end of inflation, which is determined by \M% V"(&)/V(Q)\ w 1 (since afterwards 4> 
rolls quickly to make e = 1), is given by 



-6xlO-Hv*r 1/4 (^)~ 1/2 M p 



'1/2 

$ e - $ c « -6 x 10-° (j?*) - *'* ( ) Mp, (43) 



which means that inflation ends very close to $ c (by construction) and far away from scales of interest of Eq. d42l . This places 
a constraint on the steplike function S, since the slow roll behavior should not be disrupted by the slope of this steplike function 
before the required number of efoldings of inflation. For example, if S(&) = (1 — tanh[a(<f> — <& c )])/2, a is bounded by 



M p ,<p V 100 

It is also straightforward to show that in order for $ kination domination to occur at the end of inflation, the potential that kicks 
<f> must satisfy the following condition: 

- ~ 6 ' (45) 

a much less stringent constraint than Eq. (144-b . Finally, as we have commented near Eq. d33l , we now see explicitly that if we 
make S < 10 -7 for |$ C |/M P > 10, 8 does not change the inflationary analysis. 

We have checked the details of this inflationary scenario with explicit numerical computations. The analytic discussion above 
is in good agreement with the numerical results. 



A. Inflation to Kination Domination 

The steplike feature in the potential presented in Eq. ( l33T > suggests an interesting new classical equation of motion which is 
exactly solvable that is relevant for this class of scenarios. If gravity is turned off and the In term is neglected, the equation of 
motion for the coherent state of $ at the end of inflation can be written as 

$ _ -^sech 2 (a($ - $ c )) = 0, (46) 

in which we have omitted the negligible contribution proportional to QaPc- This class of potentials allows us to now rewrite this 
equation in the limit |a<i> c | 3> 1 (the limit of interest according to Eq. d44i i) as follows: 



!> - V Q 5{<P - $ c ) w 0. 



(47) 
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The parameter a has disappeared from the equation of motion; this is not surprising given the stepfunction-like behavior of the 
potential in this limit. The beauty of Eq. ( |47] i is that it can be solved exactly. The solution is given by 

$ = e(t c - + (t - + e(t - tc){$ c + (t - t c ) 7 }, (48) 

where ©(&) is a stepfunction, 7 = ^/ 2Vo + <i>i and $j are initial values of the field and its velocity at time ti, and £ c is the 

time at which $ reaches $ c . Therefore, since the initial time variation of the field $j is small compared to \JVq , the field obtains 
a strong kick at t = t c to obtain the final state velocity of $ ~ \Z2Vq. 



B. Does Dark Matter Ever Reach Equilibrium? 

Another difference from more traditional inflationary scenarios is that the dark matter is initially out of thermal equilibrium af- 
ter reheating (or more accurately, entropy production) in this scenario. To see this, first consider the more traditional inflationary 
paradigm. The neutralino \ self-annihilation rate behaves as 

T x ~ (*v)n% ~ a 2 T, (49) 

where \ is relativistic, and a = g^ v /(4n) is the weak coupling expansion parameter. On the other hand, the expansion rate after 
reheating in standard inflationary scenarios behaves as 

H ~^W,> (50) 

which means that the neutralinos are in equilibrium as long as 

T<2.6x 10 14 Gevf^Vf— r 1/2 - (51) 

Vi/3oy uoo/ 

In our scenario, however, H is governed by $ (not T) during kination domination, which is also when \ freezes out: 

3 



H 



p , ,, ,. ., , 



In the above, we have used the fact that $ 2 a 1 /a 6 during kination domination. Hence, the neutralinos reach equilibrium 
through self-annihilations only for 

v \ 1 /' 2 

Vq 



T ^ 8Xl0GeV 'l730J U3.9xl0W ' • ^ 



Comparing this temperature with the original relativistic species temperature of 

6 x 10 7 GeV, (54) 



H e ( V ^ 1/2 



2tt V[3-9 x 10 13 GeV] 

we see that there is a long period in which the neutralino self-annihilations are out of equilibrium after the entropy in the universe 
has been produced. The fact that Eq. d53l > exists for T> 1 GeV is important, since our goal is to embed the modified dark 
matter freeze-out scenario in which the dark matter was originally in equilibrium. 



IV. OTHER PREDICTIONS FOR COSMOLOGY 



One of the main advantages of embedding the kination domination scenario within an inflationary cosmological setting is that 
other predictions for observables can be made whose experimental confirmation would either support or rule out the scenario. 
We have already seen that if tensor perturbations are measured in the near future, this scenario is ruled out. There are many other 
observables correlated with this dark matter scenario. Although an exhaustive analysis of these signatures is beyond the scope 
of this paper, we briefly discuss several possibilities here. 
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In the present scenario, an out of equilibrium effective "temperature" scale as high as lO 8 ?^ 1 ^ 4 GeV is reached at the end of 

inflation and the equilibration of the (MS)SM particles occurs by temperatures of about lO 6 ^ 1 ^ 4 GeV. During the time until 
the temperature is below ~ 1 GeV, the Hubble expansion rate differs significantly from the usual radiation dominated universe 
value. This implies that any physics which depends both on the temperature and the Hubble expansion will be modified. 

One testable example is the gravity wave production during the electroweak phase transition (see for example rf32l l33l l34ll35i 
Hit]), where the peak frequency of the gravity waves are set by the Hubble scale. Since H is about 10 5 y / 77$ of the usual 
Hubble value (for the same temperature), the effects can be large even for a very small 77$. Exact details will require a careful 
reanaiysis of the gravity wave production. The effects on electroweak baryogenesis are expected to be weaker since the out of 
equilibrium condition is primarily provided by the bubble wall velocity [38, 39]. A more careful investigation of this issue is left 
for future work. 

Another prediction of this scenario is that since the temperature at which equilibrium is reached for any heavy lepton number 
carrying particle (such as a right handed neutrino) will be relatively low as in Eq. d53l , possible leptogenesis mechanisms in 
this context will necessarily be either nonthermal or nonstandard. It will be interesting to explore what kinds of leptogenesis 
scenarios are viable for this class of models. 

There are also predictions associated with particle astrophysics. For example, it is well known that explaining the HEAT 
measurement (and other cosmic ray positron measurements) of excess positrons around and above 7 GeV requires an efficient 
annihilation of neutralinos (or other thermal relics in the context of models of extra dimensions, etc.) within our halo lf40l l4ll 
l42ll . One of the many problems associated with this efficient annihilation scenario is that the relic abundance of neutralinos 
is generically too low (by a factor of 10 to 100) to explain most of the dark matter energy within the context of standard 
thermal scenarios. The kination scenario can clearly give the necessary boost to resurrect the neutralino dark matter annihilation 
explanation of the excess positrons. There are also other cosmic ray signatures which may shed light on nonstandard H behavior 

Ha 

Finally, there are other possible signatures such as the change in BBN due to the effects of residual annihilations after freeze- 
out ll45ll . and the change in cosmic string generated gravity wave signature ll46ll47ll due to the change in T/H involved in the 
scaling behavior of cosmic strings. Since much of cosmology is about studying the out of equilibrium phenomena generated by 
the expansion of the universe within a finite temperature setting, further signatures related to a non-standard relationship between 
T and H will appear as we learn more about the early history of our universe. 



V. CONCLUSIONS 



The possibility of a new scalar is generic in extensions of physics beyond the standard model of particle physics and cosmology 
(particularly in those containing a dilatonic field degree of freedom). It is natural to expect that the dark energy density is 
connected with such a new scalar field degree of freedom. However, the elusive nature of dark energy, and its requisite small 
couplings to observable fields, make such conjectures difficult to verify or disprove. 

We have constructed a viable class of inflationary scenarios which exhibit a period of kination domination after inflation 
where the inflaton plays the role of the quintessence. Such scenarios have the intriguing feature that they can lead to observable 
consequences for the dark matter freeze-out of thermal relics expected in many TeV-scale extensions of the SM, which can be 
tested at the LHC and ILC. We have focused here on supersymmetric scenarios as prototype ex amp les, for which the connections 
between astroparticle and collider physics have been extensively explored |49[ |5oj IH [52j [53j IM[55|j56l [57 j 1^ I5§j 6 pj 6ll |62j 
l63i l64l l65l l6ol . However, other WIMP candidates have recently emerged in models with flat 1 67 , 68 ] or warped B69ll70ll extra 
dimensions, in Little Higgs theories fnl p72l p73l p74l p75ll . or in technicolor models ll76ll . to which the cosmological scenarios 
described in this paper could also be applied. 

The advantage of embedding kination-dominated quintessence models within an inflationary context is that it allows for 
many other correlated cosmological predictions which can corroborate or rule out such kination dominated scenarios. The most 
robust, nearly model independent, signature of this class of models is the absence of measurable tensor perturbations, such that 
any positive detection of tensor perturbations in upcoming experiments can rule out this scenario (at least as far generating an 
observable shift in the DM abundance is concerned). Other examples are measurements of gravity waves from the electroweak 
phase transition, shifts in the predictions for baryogenesis/leptogenesis, implications for the cosmic ray flux from dark matter 
annihilations, and other phenomena which depend on the ratio of the photon temperature to the Hubble expansion rate T/H. 

The class of models in this paper is meant to be illustrative and represent early attempts at model-building, and as such certain 
features are not optimal. For example, one might argue that the toy model presented here may not be easily achievable from an 
effective field theory point of view. However, it is reasonable to believe that this class of models are as potentially viable as any 
quintessence and most inflationary models that are considered seriously in the current literature. Given that such scenarios may 
be testable through their potentially dramatic interconnection with dark matter predictions and TeV-scale particle physics, they 
represent an intriguing and potentially fruitful ground for quintessence model-building which warrants further exploration. 
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Appendix A: PARTICLE PRODUCTION COMPUTATION 



In this Appendix, we consider the particle production computation for a real scalar degree of freedom coupled to gravity. We 
start with the perturbative expression for the Bogoliubov coefficient as an integral over conformal time [48]: 
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(A2) 



which assumes that the effective mass is dominated by the Ricci scalar R ^> m 2 , (note that £ = corresponds to minimal 
gravitational coupling). Using the approximation 



H 2 = 



V /(3M 2 ), v <7 lc 



we find 
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in which 9(z) is a unit stepfunction which evaluates to 1 for z > 0. For the limits of the integral over w 2 in Eq. ( lAlb . we have 



during inflation, and 



(l+0a 2 e ^ + k 2 < W 2 <k 2 



(A6) 



(A7) 



after the end of inflation. The conformal time r] is not a single valued function of w 2 at the transition at the end of inflation. 
However, since the non-single valued time period is arbitrarily short (in the limit that the potential behaves like a step function), 
this time period can be excised from the computation without loss of numerical accuracy, as long as a UV cutoff is imposed. 
The reason for the UV cutoff is that the peak strength of the nonadiabaticity responsible for particle production is precisely 
determined by the detailed gravitational dynamics of the transition time period, which we excise to simplify the computation. 
Since the particle production is through gravitational curvature, there will generically be an exponential cutoff in momentum of 
the particles produced at H e , the expansion rate at the end of inflation. With this simplification, Eq. ( IAU can be written as 



At 



1 f Xl dx 
- — exp 

4 J k 2 x 



-II J 2Jx~Vk 2 - x „ 
-i2V3\ - +£ I v 2arctan 



G 



dx 



exp 



x — k 2 



Vk 2 



k 2 



C^k 2 



2 In 



\/ x — k 2 + 



C 2 (k 2 



, \/ x — k 2 

in which C12 are constant phase factors which are independent of the integration variable x and depend on k 2 , while 
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which correspond to w 2 just before and after the end of inflation (approximately a step function transition). 

To obtain an estimate for \(3k | 2 , we will neglect the interference term in \0k | 2 , which means that we can neglect the integration 
constants C;. We will also utilize the fact that the contribution to the dx integral in Eq.( lA8b is appreciable only when x is near 
xi and far away from k 2 because of the damping phase oscillations near x « k 2 . Finally, we will only account for contributions 
with w 2 > 0, so as to maintain the particle interpretation of the massless modes produced. We then find 



in 



M 2 ^6 V fc 2 



n \ 2 



hi 



2% (1 W 
M? V 6 / k 2 



Q(k-k min )Q(aH e X-k), (All) 



where k min = w 4Vb / (M 2 )(l/ 6 + £) is an appropriate infrared cutoff imposed to maintain the particle production interpretation 

of the massless modes, and aH e \ is the UV cutoff, with A ~ O(l) reflecting the uncertainty in the UV cutoff function. 

Taking the approximation that the log factors contribute 0(1), the energy density of the real scalar degree of freedom that is 
produced is 
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Taking the parameterization of the thermal equilibrium abundance of one real scalar degree of freedom as p — (7r 2 /30)T c 4 ff , 
where T c s is the effective temperature and counts the number of degrees of freedom, we find 



T cff ~§(2A). 



(A13) 



For an order of magnitude estimate, we will absorb the uncertainty factor 2A in the uncertainty in the effective number of degrees 
of freedom during the "reheating" stage at the end of inflation, and therefore write 



He 

2ir 



(A14) 



throughout the paper. 



